Comments on a bound state model for a two body system 
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We show that in classical mechanics, as well as in nonrelativistic quantum mechanics the equation 
of the relative motion for a two-body bound system at rest can be replaced by individual dynamical 
equations of the same kind as the first one, but with different parameters. We assume that in rela- 
tivistic quantum mechanics the individual equations are Dirac equations with modified parameters 
in agreement with the individual Schrodinger equations. We find that products of solutions to the 
individual equations with correlated arguments are the quantum analogues of the classical repre- 

, sentation of a bound system and represent suitable models for the bound state wave functions. As 

validity test for the new representation of bound states we suggest to use some observable differences 

' between this one and the representation in terms of relative coordinates. 
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I. INTRODUCTION 

4 i In the textbooks for classical mechanics and nonrelativistic ouantum mechanics (see for instance Q and 0) the 
approaches to the two body problem start with the replacement of the dynamical equation by two equations: the first 
t-H ' one describes the free motion of the center of mass while the second equation describes a particle with the mass equal 
to the reduced mass, moving in an external potential well. A long experience has shown that the solution of the last 
one offers a suitbale representation of the bound state in the case of the heavy-light systems at rest, where it may be 
assumed that the position of the heavy particle coincides with that of the center of mass and the relative motion may 
be seen like the motion of the light particle with respect to the center of mass. 

In all the other cases the relative position vector is not sufficient for an adequate representation of a bound system. 
Moreover, the separation of the center of mass from the relative coordinates leaves the intrinsic properties of the 
particles, like spin or charges, outside the dynamical scheme, because these ones cannot be associated to either of 
them. Therefore, whenever the intrinsic properties play a significant role, the bound system must be represented in 
terms of the individual coordinates with respect to the center of mass which is supposed fixed. This is particularly 
rj - * important when the internal particles are individualized by their specific behaviour with respect to an external force. 
This is for instance the case of semileptonic decays of the heavy-light mesons where only the heavy quark is decaying, 
the other being spectator. 

In this paper we present a procedure of introducing the individual coordinates and equations of motion in the case 
of a bound system made of two bodies with arbitrary masses. We start from the observation that in the rest frame of 
the bound system the relative position vector and the individual position vectors of the internal bodies with respect 
to the center of mass are collinear. Using this fact we transform the equation of the relative motion into individual 
equations of the same kind as the first one but with different parameters and express the solution to the bound state 
problem in terms of individual solutions. 

In the next section the procedure is applied to nonrelativistic classical systems. We recall the well known fact that 
the representation of the bound system in terms of the relative coordinates is equivalent with the representation in 
terms of correlated individual coordinates. 

In the third section the procedure is extended to nonrelativistic quantum systems. From the Schrodinger equation 
for the relative motion we derive similar individual equations with modified parameters and write the complete wave 
function of the bound system as a product of individual wave functions and of the localized solution of the free wave 
equation describing the motion of the center of mass. The procedure gives rise to some specific effects which can be 
used as tests of validity. 

The case of relativistic quantum systems is discussed in the fourth section. Recalling that relativistic covariance 
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is incompatible with the separation of the relative and individual coordinates we consider that the Dirac equation 
with interaction potential must be seen as an individual equation. The parameters are chosen in agreement with the 
individual Schrodinger equation of the kind defined in the previous section. As a test of validity we calculate the 
mass and the energy levels of the hydrogen atom and notice the coincidence with the values obtained with the usual 
procedure up to terms of order or. 

In the last section we comment briefly on the relation between the coordinate and momentum representations of a 
bound system. 

Throughout the paper we use the following notations: r\, r*2 and r = r\ — r*2 are the individual and relative position 
vectors in an arbitrary coordinate system. The center of mass position vector is defined as: 

tcm = 771/1 + 772F2 (1) 

where 771.2 = mi+m 2 w ^h mi, 2 the masses of the internal particles. Obviously, 

± 7=1,2 = r C M + T]2,lf (2) 

and hence, whenever the center of mass coincides with the origin of the coordinate system one has 

r= — = . (3) 



II. CLASSICAL MECHANICS 



In classical mechanics the Newton equation of motion for the particle with reduced mass is: 

d 2 r(t) 

nmimi + mz)—^- = -V r V(r(t)) (4) 

where V(r) represents the interaction potential. 

Using the relations JQ) and © which are valid in the center of mass frame in the case of an isolated system even 
when r, f\, f 2 are functions of t, we observe that the equation Q can be written in either one of the two forms: 

mi ^W i = -^V^Ffe-VO (5) 

or 

m 2 ^^-=- V iV r2 V( V ^r 2 ). (6) 

We note that in the heavy- light case, where 772 ~ the heavy particle is at rest and the light one turns around it, as 
expected. 

The solutions to J5J and JSJl describe the individual motion of the physical particles with the masses mi and m.2 in 
the external potentials 772 F (772" ^i) and r]iV(r]i r%) which are of the same kind as V(r). If fcM — one can say that 
the instantaneous positions are at the ends of a segment of variable length, turning around the center of mass which 
is fixed and divides the segment in two parts in the ratio —r\2"f]\ X ■ Noticing the similarity of equations of motion (J3J) 
and © we also can say that the trajectories are similar and the instantaneous momenta of the two particles are equal 
and opposite. 

Now resorting to © it can be easily shown that the sum of the total energies of the individual particles is equal 
with the energy of the nonphysical particle with reduced mass. One finds indeed: 



Ei = ~^[a) +^(^iW)^2£ (7) 
E2 = ^(^) 2 + mV{ ^ r2m=rilE (8) 



where 

E = 2 + V (r(t)). (9) 

Similar relations hold for the angular momenta and hence one may conclude that in classical physics a bound system at 
rest may be represented either in terms of the relative coordinates or of the individual coordinates with the constraint 

©■ 
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III. NONRELATIVISTIC QUANTUM MECHANICS 



In quantum mechanics the internal state of a bound system is supposed to be described by a solution with finite 
norm, ip(r), of the time independent Schrodinger equation having the generic form: 

(H - E)ip(r) = (10) 

where 

H{r)=-- -r- rV' + Vff). (11) 

E is the eigenvalue of Ti. corresponding to tp an d represents the quantized value of the internal energy of the bound 
system which is the sum of the kinetic and potential energies of a particle with reduced mass in the external field 
V(r). 

Like in classical mechanics, by using @, the dynamical equation for the relative motion l|l(J|) of spinless particles 
can be written in terms of the individual coordinates either as: 



(-^- V 2 ri + 772 V^n) - m E] Mn) = 
\ 2mi / 



(12) 



or as 



7 ll , „ ^/„-i 



2TO2 



V^ 2 + m V{^ L r 2 ) - m E)Mr2 = 0. (13) 



Obviously, the equations (|1(J|) , i|12|) and 113|) with the same kind of boundary conditions are nothing else but different 
ways of writing the same Schrodinger equation. However, it is worth noticing that the solutions to the last two 
equations, tpm and "0n 2 , represent two physical particles with masses mi and m 2 individually bound in the attractive 
external potential wells 772 F^^Vi) and 771 ^(77^ 1 r 2 ) respectively. (We remark that their distance to the corresponding 
centers of forces are r\ 2 x r\ = 77^ 1 r 2 = r.) In view of the agreement with the classical case where the individual orbits 
are similar, we assume that the wave functions ip(r), , ril (r 1 ) and ipn^if^) are similar functions of coordinates which 
means that their spatial quantum numbers are equal. Also, as it results from (|12fl and l|13|) . the corresponding energy 
eigenvalues are E, E x = 772-^ and E 2 — r\\E. 

We notice that the difference between tp, ip ni , and ip„ 2 comes from the different values of the parameters in the 
interaction potential and from the different units of length used in the r, ri, and f 2 spaces. Like in classical case the 
units of length satisfy the relation: u r : u ri : u r2 = 1 : 772 : 971 . 

Now taking into account the individual character of the wave functions ip ni and tp n2 we assert that the wave function 
of the bound system at rest can be written as follows: 

f 2 ) = ipm (n - r C M)^n 2 (r 2 - ?cm) A a (r C M) (14) 

where the function A CT is the solution of the Schrodinger equation for the center of mass which is localized in a region 
of range tr _1 around fcta = 0. In the limit a — > 00, „ v ^ is the probability density to find the two particles at the 
ends of a segment of length \(f\ —v%)\ divided by the center of mass in two parts in the ratio — 

We remark that 1 / 1 , I / 2 (ri, 7*2) satisfies the main conditions to be considered a suitable model for the bound state 
wave function: its L 2 norm in the B? x i? 3 space is finite and in the limit a — > 00 the representation of the bound 
system in the coordinate space is in agreement with the classical representation. 

However, we notice that, unlike the classical case, where the correlation of the individual positions in the rest 
frame entails the correlation of the individual momenta, in the quantum case the perfect correlation of the individual 
positions transforms into a total lack of correlation of the individual momenta. Indeed, if we take A CT (rc , j\/) = 

^-^=^ exp[—r 2 - ;M <j' 2 ] and then expand the individual wave functions in terms of momentum eigenstates and integrate 
over the degrees of freedom of the center of mass we find 

^rcM^in,^) = (2ir) 3 /d 3 fc 1 d 3 fc 2 d 3 Q7/;„ 1 (fc 1 )e- l£lFl 7/;„ 2 (fc 2 )e-^ F2 A(Q) ^(fcj +k 2 -Q) (15) 



where ip(k) = (2n) 3 J d 3 rip(r)e lkr and A = (^^/2~Ka) 2 exp[— -^{pi + p 2 ) 2 }- It follows that the perfect correlation 
of the individual positions (tcm = 0) is achieved for a — > 00, while the perfect correlation of the individual momenta 
(fci + k 2 = 0) results for a — > 0. 
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Then in order to have a certain degree of correlation both of the individual positions and of the individual momenta 
the parameter a must be finite. This means that the center of mass of the two particles is distributed in a region of 
range c -1 around the origin or, in other terms, it is not fixed when the bound system is at rest. This might appear as 
a misleading view but it could be accepted if we resort to the quantum field representation of a bound system where 
the binding forces are generated by the continuous exchange of quanta between the constituents. Then the motion 
of the center of mass may be seen as the consequence of the imperfect cancellation of momenta during the quantum 
fluctuations which generate the binding. 

Closing this section we notice that in the new representation of a bound system there are two individual wave 
functions instead of a single function of the relative coordinates. Accordingly, a term in the relative wave function 
behaving like r l when r — ► is replaced by r 21 in the limit a — ► oo. Also, if the angular momentum had the eigen 
value I in the old representation, now it may have any value in the range [0, 21]. Noticing the direct influence of these 
changes on the form factors we consider that a comparison of the calculated values with the experimental data may 
provide a good test of validity for the present model. 



IV. RELATIVISTIC QUANTUM MECHANICS 

As it is well known the separation of the center of mass coordinates from the relative ones cannot be performed in 
the case of the Dirac equation and hence we are prevented from applying the procedure developed in the preceding 
section to fermion systems. 

In exchange, we are allowed to assume that the Dirac equation with an interaction potential is one of the individual 
equations if it is in agreement with one of the Schrodinger equations written above. Then, taking i|12|) as a model, 
the suitable form of the individual Dirac equation for the particle 1 is the following: 



(-i77 2 V ri <3i + /3iTOi + 772^(772 n.) - Ei)ip ni (fi) = 0. 



(16) 



The second individual equation may be obtained from (|16[) by interchanging the indices 1 and 2. 

The wave function of the bound fermion system is supposed to be of the same form like l|14[) where this time the 
individual wave functions ip ni 2 are bi-spinors and A is simply a function constraining the center of mass degrees of 
freedom in agreement with the nonrelativistic case. The mass of the bound system is the total energy of the bound 
system in its rest frame and hence it is the sum E\ + E 2 of the eigenvalues of the Dirac Hamiltonian which include 
the individual rest masses. 

A simple test of validity of the present model is to compare the difference T> between the energy levels of a hydrogen- 
like atom using the new and the old prescriptions. Specifically, one has to calculate the difference between the energy 
of the bound system at rest written as the sum of the individual energies of the electron and of the nucleus in Coulomb 
potentials with modified strengths as it results from (|16|) and the energy obtained by solving the Dirac equation for 
a particle with reduced mass in the Coulomb potential. 

Considering the expressions quoted in Ref.|3j] for the energy levels, in the first case we have 



M„. 



E\ + Ei 



(17) 



where 



E lj2 = m 



1.2 



(772, iZa) 2 



(n 



(18) 



with 



°1,2 



I-2 



3 + 



- (772, iZa) 2 



(19) 



The energy of the particle with reduced mass in the Coulomb field of the nucleus is obtained from E\ by replacing 



'12 



1 and 



777,1 



(mi+m 2 ) 

the energy levels becomes: 



. Then, after including the recoil corrections (see the difference between the values of 



V 



Z a {mi + 7772) 



8??. 



1 - 



777 1 7772 
(77li + 777 2 ) 2 



777,1 "*2 



3(t77i 



m 2 ) 



2j + 1 



7771 ™2 



(mi 



m 2 ) 



(20) 



5 



where Z is the nucleus charge. The terms in the first square brackets in (|20(l represent the difference between the 
recoil corrections and can be measured in principle. We mention that in the old treatment the recoil effect has been 
obtained by using some particular assumptions, while in the present formalism it results from the individual treatment 
of the particles. The last part is the difference between the fine splitting terms which is negligible small in the case 
of hydrogen- like atoms where — << 1. 

Closing this section we notice an important consequence of the individual, relativistic treatment of the internal 
particles: the attractive forces create a positive mass defect 5m = (mi + ^2) — M a which tends to when Za — > 
and is automatically included into the mass of the bound system together with the free masses. From the formal 
point of view this is perhaps the most impressive result of the present model. 

V. COMMENTS AND CONCLUSIONS 

In the last two sections we have shown that the procedure leading to a representation of a bound system in terms of 
individual coordinates outlined in classical mechanics also may be extended to quantum physics with similar results. 

However, it is worth noticing some of the pecularities of the quantum representations. A first one is that in 
quantum mechanics the perfect correlation of the individual coordinates is incompatible with the perfect correlation 
of the individual momenta and conversely. This is unlike the classical case where the correlation of the individual 
position entails the correlation of the individual momenta. 

A second pecularity concerns the dependence on time of the bound state wave function in the coordinate and in the 
momentum representation. In the first case the wave function is stationary and its dependence on time concentrates 
in the overall factor exp[i(Ei + E2)t]. In the momentum space the bound system is represented by the projections 
of the wave function on two particle free states whose dependence on time is exp[i(ei + e.2)t] where e% t 2 are the 
relativistic free energies. Observing that in the relativistic case E\ + E2 = M, while e\ + e2 can take any value 
in the range [(mi + 1712), 00] one may conclude immediately that a two body bound system cannot be adequately 
represented in terms of two body free states. This is not unexpected, because neither the potential energy of the 
bound system, nor the correlation of the individual coordinates can be explained with their aid. In principle this 
ones ought to be explained in field theory but, if one wishes to obtain a suitable representation of a bound system in 
the framework of the quantum mechanics, it is convenient to accept the existence of a neutral, vacuum-like, effective 
constituent beside the " valence" particles [4| , whose features compensate the differences between the momentum and 
the coordinate representations. Specifically, the momentum of this one must compensate the free individual momenta 
and its contribution to the energy must supply the system with the potential energy which misses in the momentum 
representation. 

Obviously, the effective constituent is not an elementary excitation of the background field and it is hard to 
imagine its representation in the coordinate space. However, taking into account its role in the bound system, it may 
be associated with the external potential which generates the forces binding the valence particles together. Then, 
recalling that in field theories binding is the result of the continuous exchange of quanta, we consider reasonable to 
assume that the effective constituent is an average representation of countless quantum fluctuations. 
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